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ABSTRACT 

The process of cosmic ray acceleration up to energies in excess of 10^° eV at relativistic 
shock waves with large Lorentz factors, F ^ 1 requires ~ particle energy gains 
at single reflections from the shock (cf. Gallant & Achterberg 1999). In the present 
comment, applying numerical simulations we address an efficiency problem arising 
for such models. The actual efficiency of the acceleration process is expected to be 
substantially lower than the estimates of previous authors. 

Key virords: acceleration of particles - UHE cosmic rays - gamma ray bursts - shock 
waves 



1 INTRODUCTION 

In attempt to substitute a single question mark for the 
previous two, some authors try to identify the process ac- 
celerating particles to ultra high energies ('UHE', energy 
E > 10^** eV) with ultrarelativistic shock waves considered 
to be sources of 7-ray bursts (cf. Waxman 1995a,b, Vietri 
1995, Milgrom & Usov 1995, Gallant & Achterberg 1999, = 
GA99). In the proposed models, reaching cosmic ray energies 
m excess of 10^" eV requires that particles are reflected from 
the shock wave characterized with a large Lorentz gamma 
factor, r, to enable a relative energy gain - in a single re- 
flection - comparable to (cf. GA99). On the other hand, 
basing on our experience with numerical modeling, we sug- 
gested (Bednarz & Ostrowski 1998, Ostrowski 1999) that 
such processes cannot actually work due to low efficiency of 
particle reflections. 

In the present note we elaborate this problem in detail 
with the use of numerical modeling of particle interactions 
with the shock. We show that T^' reflections can occur with 
a non-negligible rate only if a substantial amount of turbu- 
lence is present downstream of the shock. However, even in 
such conditions the number of accelerated particles is a small 
fraction of all cosmic rays hitting the shock. 



2 SIMULATIONS 

As discussed by Bednarz & Ostrowski (1998), and in detail 
by GA99, particles accelerated in multiple interactions with 
an ultrarelativistic shock wave gain on average in a single 
'loop' - upstream-downstream-upstream - the amount of 
energy comparable to the original energy, < AE >« E. It 
is due to extreme particle anisotropy occurring in large P 
shocks. Particles hitting the shock wave the first time, with 



their isotropic upstream distribution, can receive higher en- 
ergy gains. In this case an individual reflection from the 
shock may increase particle energies on a factor of ~ P^. 
However the effectiveness of such acceleration depends on 
how many of particles from the original upstream popula- 
tion can be reflected. To verify it we use Monte Carlo sim- 
ulations similar to the applied earlier to derivation of the 
accelerated particle spectra (Bednarz & Ostrowski 1998, for 
details see also Bednarz & Ostrowski 1996). The code re- 
produces perturbations of particle trajectories due to MHD 
turbulence by applying discrete scatterings of particle direc- 
tion within a narrow cone along its momentum vector. A 
procedure uses a hybrid approach involving very small scat- 
tering angles close to the shock and larger angles further 
away from it. Between successive scatterings particle trajec- 
tories are derived in the uniform background magnetic fleld. 
The respective scaling of the time between the successive 
scattering acts close and far from the shock mimics the same 
turbulence amplitude everywhere. The scattering amplitude 
was selected in a way to reproduce a pitch angle diffusion 
process for particle momentum. It requires the angular scat- 
tering amplitude of particle momentum vector to be much 
smaller than the particle anisotropy. All computations were 
done in a respective local plasma rest frame and the Lorentz 
transformation was applied to every particle crossing of the 
shock. 

Particles with initial momenta po taken as a momentum 
unit, Po = 1, were injected at the distance of 2rg (r-g - par- 
ticle gyroradius) upstream of the shock front. For all such 
particles we derived their trajectories until they crossed the 
shock upstream, or were caught in the downstream plasma, 
reaching a distance of 4rg downstream of the shock. For 
each single particle interaction with the shock the particle 
momentum vector was recorded, so we were able to consider 
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angular and energy distributions of such particles. We con- 
sidered shocks with Lorentz factors F = 10, 160 and 320. For 
each shock we discussed the acceleration processes in con- 
ditions with the magnetic field inclination ^ = 0° , 10°, 70° 
and with 16 values for the turbulence amplitude measured 
by the ratio t of the cross-field diffusion coefficient, k±, to 
the parallel diffusion coefficient, /ty. The applied values of 
r were taken from the range of (3.2 • 10~®, 0.95), approx- 
imately uniformly distributed in logT . In each simulation 
run we derived trajectories of 5-10'' particles with initial mo- 
menta isotropically distributed in the upstream rest frame. 



3 EFFICIENCY OF T^' REFLECTIONS 

In the downstream plasma rest frame the ultrarelativistic 
shock moves with velocity c/3. This velocity is comparable 
to the particle velocity c. Therefore, from all particles cross- 
ing the shock downstream only the ones with particular mo- 
mentum orientations will interact with the shock again, the 
remaining particles will be caught in the downstream plasma 
fiow and advected far from the shock front. In the simular 
tions we considered this process quantitatively. However, let 
us first present a simple illustration. 

Large compression ratios occurring in ultrarelativistic 
shocks, as measured between the upstream and downstream 
plasma rest frames, lead for nearly all oblique upstream mag- 
netic field configurations to the quasi-perpendicular config- 
urations downstream of the shock. Thus, let us consider for 
this illustrative example a shock with a non-perturbed per- 
pendicular downstream magnetic field distribution. Particle 
crossing the shock downstream with inclination to the mag- 
netic field 9 and the phase (j> - both measured in the down- 
stream plasma rest frame, = 7r/2 for particles normal to 
the shock and directed downstream - will be able to cross 
the shock upstream only if the equation 



I' 



: rg [cos(^ -|- Wgt) — cos 0] 



(3.1) 



j§ sin 6 is the 



has a solution at positive time t. Here r. 
particle gyroradius, cjg = ^ is the gyration frequency, and 
other symbols have the usual meaning. 

An angular range in the space {9, 4>) enabling particles 
crossing the shock downstream to reach the shock again can 
be characterized for illustration by three values of 9. Par- 
ticles with sin^ = 1 are able to reach the shock again if 
€(1.96, 3.48), with sm9 = 0.5 if g(2.96, 3.87) and with 
sin^ = 1/3 only for = 4.71. That means that all particles 
with (j) smaller than 1.96 (Fig. 1) are not able to reach the 
shock again if fluctuations of the magnetic field downstream 
of the shock are not present. 

For perturbed magnetic fields some downstream trajec- 
tories starting in the {9, (p) plane outside the reflection range 
can be scattered toward the shock to cross it upstream. We 
prove it by simulations presented in Fig. 2. One may ob- 
serve that increasing the perturbation amplitude leads to 
increased number of reflected particles, reaching « 13% in 
the limit of t = 1. For large magnetic field fluctuations the 
mean relative energy gains of reflected particles are close 
to 1.2r'^ for the shock Lorentz factors considered. One may 
note that for small ip and F the energy gain increases with 
growing r. The points resulting from simulations for the 




Figure 1 . A distribution of particle phases for particles crossing 

the shock downstream (as measured in the downstream plasma 
rest frame), if their upstream distribution was isotropic. A dashed 
line delimits a range of particle phases below which particles are 
not able to reach the shock again at the perpendicular uniform 
downstream magnetic field. 
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Figure 2. A ratio of the number of reflected particles, n, to all 
particles crossing the shock downstream, no, as a function of the 
magnetic field fiuctuations amplitude, t. 
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Figure 3. Momentum gains of reflected particles, p/po, as 
a function of the magnetic flcld fluctuations amplitude, r. For 
large magnetic field fluctuations tiic momentum gain approaches 
1.2r^ independently of the shock Lorentz factor. 
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smallest values of r were not included into Fig. 3 because of 
small number of reflected particles (cf. Fig. 2). 



4 SUMMARY 

We have shown that efficiency of T^' reflections in ultra- 
relativistic shock waves strongly depends on fluctuations of 
magnetic field downstream of the shock. In the most favor- 
able conditions with high amplitude turbulence downstream 
the shock the reflection efficiency is a factor of f or more 
smaller than the values assumed by other authors. More- 
over, due to the magnetic field compression at the shock, 
we do not expect the required large values of k±/k\\ to oc- 
cur behind the shock (cf. a different approach of Medvedev 
& Loeb 1999). Therefore, with the actual efficiency of 1 - 
10 % there is an additional difficulty for models postulating 
UlfE particle acceleration at GRB shocks (cf. GA99). Let us 
note, however, that the mean downstream trajectory of the 
reflected particle involves only a fraction of its gyroperiod. 
Thus the presence of compressive long waves in this region, 
leading to non-random trajectory perturbations could mod- 
ify our estimates. 
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